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^sO , New observations with atmospheric neutrinos from the underground experiments 

SuperKamiokande, Soudan 2, and MACRO, together with earlier results from 
£NJ ' Kamiokande and 1MB, are reviewed. The most recent observations reconfirm aspects 

of atmospheric flavor content and of zenith angle distributions which appear anomalous 
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■ in the context of null oscillations. The anomalous trends, exhibited with high statistics 

in both sub-GeV and multi-GeV data of the SuperKamiokande water Cherenkov ex- 
periment, occur also in event samples recorded by the tracking calorimeters. The data 
are well-described by disappearance of flavor neutrinos arising in oscillations with 
dominant two-state mixing, for which there exists a region in (sin 2 29, Am 2 ) allowed 
by all experiments. In a new analysis by SuperKamiokande, — > v r is favored over 
Vn — > Vsteriie as the dominant oscillation based upon absence of oscillation suppression 
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D , from matter effects at high energies. The possibility for sub-dominant za, — ► v e oscilla- 

tions in atmospheric neutrinos which arises with three-flavor mixing, is reviewed, and 
intriguing possibilities for amplification of this oscillation by terrestrial matter-induced 
■ resonances are discussed. Developments and future measurements which will enhance 

^ , our knowledge of the atmospheric neutrino fluxes are briefly noted. 

1 Atmospheric Neutrino Beamline 

We are lucky, you and I, to be born here on planet Earth and to have as our birthright 
the unrestricted use of a splendid neutrino beamline. Truly remarkable is that the 
originating hadronic primary beam, namely the cosmic ray flux of protons and as- 
sorted stable nuclei, is isotropic to high degree. Moreover the beamline target region, 
which is the terrestrial atmosphere, is very nearly spherically symmetric. Together 



1 Plenary talk at the XIX International Symposium on Lepton and Photon Interactions at High 
Energies, Stanford University, August 9-14, 1999. 
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these two attributes ensure that neutrino fluxes of this beamline, in the absence of 
neutrino oscillations, must be up/down symmetric with respect to the horizon |], |[. 
Consequently, observation of a sizable neutrino flux up/down asymmetry by the user 
community is clear evidence that new physics is happening with neutrinos. The null 
oscillation up/down symmetry for neutrino fluxes is however not complete. There are 
geomagnetic effects which produce mild distortions in the fluxes of low energy neu- 
trinos and of horizontal neutrinos. These distortions, which are latitude-dependent, 
provide useful tests for data verity. The beamline delivers neutrino fluxes which are 
wide-band in E v and which contain the flavors {v^+V^) and [y e +V e ). In the absence 
of oscillations these flavors must occur in the ratio 2:1 for sub-GeV neutrinos; for 
multi-GeV neutrinos the ratio should increase gradually with E v . At or below the 
Earth's surface the atmospheric flux is about 10 3 z/'s incident per human body per 
second 0] , an amount which is adequate for experimentation but does not pose a radi- 
ation safety hazard. The neutrino path lengths L which are possible in this beamline 
range from 20 km for i/'s incident from the local zenith, to 13,000 km for u's arriving 
from the opposite side of the globe. Within the beamline there are regions of dif- 
ferent, roughly uniform, matter densities. These include the Earth's mantle (density 
~ 4.5 g/cm 3 ) and the Earth's core (density ~ 11.5 g/cm 3 ). This arrangement may 
eventually permit experimental strategies to be tried which are akin to utilization 
of regeneration plates in K° beams. The neutrinos from this ever-running beamline 
give rise to useful reaction final states both in and below any detector deployed un- 
derground. By investigating the full panoply of event types possible with charged 
current (CC) or neutral current (NC) interactions, experimentalists can explore the 
physics of atmospheric neutrinos for incident E u ranging from 100-200 MeV up to 
and exceeding 1000 GeV. 

2 Oscillation Phenomenology 

We believe there to be three active neutrinos; there may be sterile ones v s as well. 
For the active neutrinos, the weak flavor eigenstates u e , u^, and v T are related to the 
mass eigenstates according to a product involving the unitary mixing matrix U : 
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The oscillation probabilities which follow from this can, in principle, involve numerous 
competing processes: 

P{u a - up) = 5 aP - E u * ■ U & ■ V U ■ U h ■ sin2 

i=l j=i+l 
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Fortunately there are cases wherein oscillations decouple so that the situation is well- 
described by two-neutrino oscillations, for which the mixing matrix is much simpler: 
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The probability for oscillation between the two participating flavors can then be 
written using the well-known expression 
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It is convenient to define the vacuum oscillation length Lq: 
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The oscillation phase can then be expressed as (ttL/Lq) 



3 The Underground Detectors 

Currently there are three underground experiments which are accumulating atmo- 
spheric neutrino data. The premier detector in this field is SuperKamiokande (Super- 
K). It is a 50 kiloton water Cherenkov detector deployed in a configuration of two 
concentric cylindrical volumes. The inner volume is the 22.5 kiloton fiducial region, 
while the surrounding outer volume is used to veto entering tracks and to tag exiting 
tracks. Flavor tagging of events is based upon the relative sharpness or diffuseness of 
Cherenkov rings, with muon tracks yielding sharp rings, and electrons yielding diffuse 
ones ||. The analyzed exposure for Super-K in-detector neutrino reactions reported 
here is from 848 livedays; this corresponds to a whopping 52 fiducial kiloton years! 

MACRO and Soudan 2 are tracking calorimeter detectors. MACRO is a large-area, 
planar tracker. It is optimized for tracking in vertical directions and is sufficiently 
massive (about 5.3 kilotons) to be effective as a neutrino detector. Charged particle 
tracking is carried out using horizontal layers of streamer tubes with wire and stereo 
strip readout. Three horizontal planes and also vertical walls of liquid scintillator 
counters provide timing information with resolution of about 0.5 nsec ||]. MACRO 
has the largest rock overburden of the three underground experiments, consequently 
the flux of downgoing muons which can give rise to backgrounds is lowest at its site 

1- 

Soudan 2 is a fine-grained iron tracking calorimeter of total mass 963 tons which 
images non-relativistic as well as relativistic charged particles produced in neutrino 
reactions. The detector operates as a slow-drift time projection chamber. Its tracking 
elements are meter-long plastic drift tubes which are placed into the corrugations of 
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steel sheets. The sheets are stacked to form a tracking lattice of honeycomb geometry. 
A stack is packaged as a calorimeter module and the detector is assembled building- 
block fashion using these modules. The calorimeter is surrounded on all sides by 
a cavern-liner active shield array of two or three layers of proportional tubes @. 
The contained event sample reported here is obtained from a 4.6 fiducial kton-year 
exposure. 

4 Atmospheric Neutrino Flavor Ratio 

A hypothesis test of long-standing for the existence of anomalous behavior of atmo- 
spheric neutrinos is "the flavor ratio" for which updated measurements have become 
available this summer. Atmospheric neutrinos are produced almost entirely in pion - 
muon decay chains initiated by cosmic ray interactions in the upper atmosphere. As 
a consequence, (i^+F^) versus (u e +V e ) neutrino flavor rates occur in a ratio 2:1. The 
underground experiments examine the ratio-of-ratios R, which is (z/ M + z7 M )/(z/ e + V e ) 
from the data, divided by the same ratio from a Monte Carlo. In the absence of new 
physics the ratio-of-ratios should be unity; and so the degree to which deviation from 
unity is observed is a measure of anomalous behavior of the fluxes. In actual practice, 
the experiments measure a related quantity, R', the ratio of observed event counts. 
For the Super-K water Cherenkov experiment, R' is the ratio of single-ring /x-like to 
e-like events in the data divided by /i-like to e-like from the Monte Carlo 0. For 
Soudan 2, R' is the ratio of single-track to single-shower events for the data, divided 
by the same ratio from the Monte Carlo ||. 

Here, then, are the latest results from Super-K, updated to include the 848 day 
exposure: For the "sub-GeV" sample (with event visible energy E vis < 1.33 GeV), 

^'(/i-ring/e-ring) = 0.68 ± 0.02 (stat.) ± 0.05 (syst.) . 
For the "multi-GeV" sample (E vis > 1.33 GeV), 

#'(//-ring/e-ring) = 0.68 ± 0.04 ± 0.08. 

From the Soudan 2 iron calorimeter there is an updated measurement based upon 
contained track and shower events of a 4.6 fiducial kiloton year exposure. The events 
occur mostly within the sub-GeV regime as defined by Super-K: 

R\trk/shwr) = 0.68 ±0.11 ± 0.06. 

Measurements of the atmospheric flavor ratio have been accumulating from the 
underground experiments for more than a decade [^, |H| [II], [12| . These most recent 
results reconfirm the atmospheric anomaly as first reported years ago by the water 
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Figure 1 : Measurements of the atmospheric neutrino flavor ratio of ratios jl3) . 

Cherenkov experiments Kamiokande and 1MB. Fig. [I] shows that the various R' 
measurements, by different experiments with different techniques and systematics, 
give a consistent picture. The flavor content of the atmospheric neutrino flux is 
anomalous but in a way that is readily understandable, if indeed muon neutrinos are 
being depleted by — > v x oscillations over pathlengths which occur in the terrestrial 
beamline. 



5 Zenith Angle Distortions and Super-K Data 

To elicit the pathlength L dependence which, in an oscillation scenario, will corre- 
late with [y^ + 77^) disappearance, we consider the distributions of neutrino zenith 
angle which have been obtained for fully contained (FC) and for partially contained 
(PC) events in the SuperKamiokande experiment. In evaluating zenith angle distri- 
butions and also flavor ratios, it is useful to keep in mind trends which are shown 
by the survival probability curves in Fig. ^|a for neutrinos [14]. The curves de- 
pict the probability for — > from an atmospheric flux for which cos #.2 at 1.0 is 
vertically downgoing and cos# 2 at -1.0 is vertically upgoing. The curves are drawn 
for "representative" v x parameter settings which we use again in paragraphs 

below, namely sin 2 29 = 1.0 and Am 2 = 5 x 10~ 3 eV 2 . The oscillation pattern in 
Fig. ^]a evolves in a regular way with increasing energy of the neutrino. For E v of 
250 MeV, the first oscillation swing severely depletes the downward-going flux, and 
rapid oscillations deplete the flux incident from below-horizon; the net result is a 
substantial average depletion at all incident angles. At energies above 1 GeV however 
the depletion moves almost entirely to the flux incident from below-horizon, and 
this situation remains for E v increasing to 30 GeV. At higher E v the pattern shifts 
to beyond range, and depletion ceases because our planet is not big enough to 
accomodate the first oscillation swing. 

Distributions showing ten bins in cos 6 Z for events of the 848-day Super-K exposure 
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Figure 2: a (left): Survival probability curves for mono-energetic, isotropic fluxes of muon neutri- 
nos with Vp, — ► v r and sin 2 29 = 1.0, Am 2 = 5 x 10~ 3 eV 2 for nine values of E v . b (right): Super-K 
distributions for cos# z single ring e-like and /i-like events. Expectations for null oscillation and for 
— > v T are shown by gray-line and solid-line histograms respectively. 

are given in Fig. |2|b. The v e and flavor samples are shown separately and are 
subdivided according to E V i S . The v e events show no angular distortion in either the 
sub-GeV or multi-GeV regimes. In striking contrast the samples show large regions 
of disappearance, the samples being depleted relative to expectations of the null 
oscillation Monte Carlo (gray-line histograms). The depletions exhibit dependence 
on zenith angle and therefore on path length L. Additionally, the depleted regions 
are of different character in the sub-GeV and multi-GeV sets. At sub-GeV energies 
the /i-like events appear depleted at all angles including those with incidence from 
above horizon. At multi-GeV energies however, the depletion is mostly restricted 
to incidence from below-horizon. Although the correlation between the final state 
lepton and the initial neutrino direction is relatively poor for sub-GeV compared to 
multi-GeV data, nevertheless the trend is suggestive of a dependence on E v for 
flavor disappearance. As shown by the solid-line histograms superposed in Fig. ^|b, 
the zenith angle distortions of the flavor samples are well-described by a fit of two- 
state Up — > v T neutrino oscillations (discussed below). Contrastingly, the v e samples 
are in agreement with the null oscillation Monte Carlo (MC) to a degree which is 
perhaps disappointing. With multi-GeV u e 's which presumably traversed the Earth's 
core, for example, no irregularity is apparent; there are no hints anywhere to suggest 
Vjx v e oscillations. 

The depletion of neutrinos can be shown in an informative way by plotting 
the asymmetry in zenith angle as a function of event momentum as in Fig. |3|a. The 
asymmetry A is defined A = (U — D)/(U + D), where U is the number of events with 
upward incidence at angles cos 9 > 0.2 and D equals events with downward incidence 
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at angles cos# < —0.2. For the single-ring e-like events, A ~ at all momenta. For 
the /i-like events however, A becomes increasingly negative, there being a dearth of 
upward-going versus downward-going neutrinos which becomes more pronounced with 
increasing momentum. For the multi-GeV sample the value of A is — 0.32±0.04±0.01, 
which is nearly eight standard deviations from zero asymmetry. 
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Figure 3: a (left): Zenith angle asymmetry versus event momentum for single ring e-like and 
/i-like events of Super-K. b (right): Super-K allowed region in the neutrino oscillation parameter 
space for — > v T mixing, based upon the x 2 fit to FC and PC single ring event distributions. 



For the fully-contained and partially contained single-ring data just shown, the 
Super-K collaboration uses a \ 2 function to determine the oscillation parameters of 
two-state mixing: 

X 2 (sin 2 29, Am 2 , e) = £ ~ + £ % (6) 

i=l,70 °i j a j 

For this purpose the /i-like and the e-like samples are sub-divided using five bins in 
cos 9 Z and seven bins in momentum. The \ 2 is the sum of data minus MC expectation 
squared over the 70 bins, where the MC is a function of the oscillation parameters 
sin 2 26 l , Am 2 , and parameters Sj which allow for systematic effects. The Sj include 
the parameter a which appears in the flux normalization factor (1 + a). At each 
point in the plane of sin 2 29 and Am 2 , the \ 2 is minimized with respect to the Ej 
parameters; the minimum \ 2 point (best fit) is then obtained. Contours for allowed 
regions at 68%, 90%, and 99% CL are obtained on the basis of \ 2 ~~ Xminimum as 
shown in Fig. ^o. The — ► v T oscillation best fit yields x 2 = 55/67 d.o.f. and fares 
much better than the null oscillation fit having 177/69 d.o.f. At the best fit point the 
oscillation parameter values are sin 2 29 = 0.99 and Am 2 = 3.1 x 10~ 3 [eV 2 ]; the MC 
flux normalization is shifted upwards (a = +0.05) relative to the absolute rate based 
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upon the Honda et al. fluxes for the Super-K site [|T5|. It is comforting to see that 
sin 2 29 from the best fit with Super-K FC and PC events now occurs in the physical 
region, for this has not always been the case in the past. 
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Figure 4: Two track-plus-recoil events, a single shower event, and a flavor multiprong event, 
recorded in Soudan 2. 



6 Contained Events in Soudan 2 

Concerning evidence for neutrino oscillations carried by in-detector neutrino interac- 
tions, a "second look" is afforded by the fully contained track, shower, and multiprong 
events recorded by Soudan 2. Projected images of "typical" data events are shown 
in Fig. |j; these include two examples of muon tracks with companion recoil protons 
[y^ quasi-elastics) , a shower event [y e quasi-elastic), and a z/^-flavor multiprong. 

The approach taken by Soudan 2 is to isolate a sub-sample of events for which 
L/E v can be measured with good resolution on an event-by-event basis, thereby 
allowing the oscillation analysis to be carried out directly using L/E v distributions. 
In a fine-grain tracking calorimeter, E u can be reliably estimated based upon E vis . To 
ensure good resolution for ascertaining the incident neutrino direction, quasi-elastic 
single track and shower events are selected which have measurable recoil protons. 
Otherwise, in the absence of a visible recoil, the lepton energy is required to exceed 
600 MeV. Multiprong events are also selected, provided that E vis > 700 MeV and 
|£p„i s | > 450 MeV/c and Pi ep t > 250 MeV/c. The momentum requirements improve 
the resolution of neutrino direction and ensure reliable flavor-tagging for charged 
current events (success probability > 0.92). For the (u e ) sample, AE/E is 20% 
(23%). For pointing of the event along the original neutrino direction, the resolutions 
are of order 20-30 degrees which is quite respectable for a sub-GeV data set ]TB| . 



Zenith angle distributions for the resolution-enhanced (HiRes) u e and v ^ samples 
are shown in Fig. |5]a, where the MC rates have been normalized to the observed 
number of v e events. For the u e sample, the zenith angle distribution follows the 
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Figure 5: a (left): Distributions in cos9 z for i/ e and flavor event samples. Data (crosses) are 
compared to the null oscillation MC (dashed histogram) where the MC has been rate normalized to 
v e data, b (middle): Distributions of \og(L/E u ) for u e and charged current events compared to 
the neutrino MC with no oscillations; the MC has been rate- normalized to the v e data, c (right): 
Comparison of L/ E u distribution for data (crosses) and expectations from neutrino oscillations 
for four Am 2 values, with sin 2 29 = 1.0. 



shape predicted by the MC without oscillations. The corresponding data distri- 
bution however consistently falls below the MC expectation. The dearth is mild but 
discernible for v incidence above the horizon and becomes more pronounced with 
below-horizon incidence. These features are in agreement with those exhibited with 
much higher statistical weight by the sub-GeV FC single-ring events of Super-K. 

Distributions in log LjE v for the HiRes v e and samples are shown in Fig. [5|b 
wherein the data (crosses) are compared to the null oscillation MC. The peaks at 
low log L/E u are populated by down-going neutrinos incident from above-horizon; 
the lower flux central regions are populated by horizontal neutrinos, while the peaks 
at higher log L/E u contain neutrinos traveling upward through the Earth. To within 
statistical fluctuations, the v e sample follows the null oscillation MC expectation. 
For the sample, there is a depletion which pervades the entire up-going region 
and extends into the down-going flux, subsiding only in the lowest L/E u bins which 
contain the most vertically down-going events. 

The log L/E v distributions from data can be fitted to oscillation- weighted MC 
events using a x 2 constructed similarly to the function utilized by Super-K. An ex- 
ploratory matchup is shown in Fig. |^c, for which sin 2 29 is set to 1.0 and the data 
is displayed together with weighted MC distributions for — > v x oscillations at four 
different Am 2 values. At Am 2 = 10 -4 eV 2 the oscillation solution exceeds the data 
in the down-going hemisphere. At Am 2 = 10~ 3 eV 2 the matchup improves, and at 
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Figure 6: a (left): x 2 "map" over the Am 2 , sin 2 29 plane; the MC normalization is also allowed 
to adjust. Contour boundaries are approximately 68% CL, 90% CL, and 95% CL based upon Xmin 
plus 2.30, 4.61, and 5.99 respectively, b (right): Soudan 2 allowed regions for z/„ — > v x oscillations 
at approximately 68% and 90% CL. 



Am 2 = 8 x eV 2 the oscillation solution follows the data rather well. However 
Am 2 = 10" 1 eV 2 , is "too far" - the oscillation solution falls below the data in the 
down-going hemisphere and in the up-going hemisphere as well. This sequence illus- 
trates key features of the x 2 mapping of the (Am 2 , sin 2 29) plane shown in Fig. |6|a. 
The best fit lies in the darkened basin region of the contour. The boundaries of the 
different shaded areas correspond to regions allowed at approximately 68%, 90%, and 
95% CL. 

In the contour map there appears a 'ridge of improbability' at Am 2 lying just 
above the best fit 'basin'. This region corresponds to oscillation solutions for which 
the first oscillation swing should create a depletion in the downward-going flux. 
No such depletion occurs in the data and since the events have sufficient directional 
resolution to show it if it occurred, the x 2 becomes large there. The projection of 
the x 2 contours onto the (Am 2 , sin 2 29) plane is shown in Fig. |6|b. The minimum x 2 
point is at sin 2 29 = 0.95 and Am 2 = 0.8 x 10~ 2 eV 2 . The flux normalization, which 
is allowed to vary in the fitting, is reset at the minimum x 2 point to 0.82 times the 
absolute event rate based upon the Monte Carlo. The Monte Carlo utilizes the 1989 



Bartol flux calculation for the Soudan site 117 



7 Partially Contained Events in MACRO 

Although the Soudan results are in general agreement with the neutrino oscillation 
effects reported by Super-K, they do not at present confirm the striking depletion in 
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Figure 7: a (left): Cross section sketch of the MACRO detector, illustrating the event topologies 
of partially contained and through-going samples, b (right): Distributions of parent neutrino 
energies for the three neutrino samples of MACRO. 



upgoing muon neutrinos shown by Super-K's multi-GeV events. Fortunately, event 
samples which provide another, independent viewing of the multi-GeV E v regime are 
being accumulated in the MACRO experiment |[L8||. 
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Figure 8: Distributions in cos8 z for MACRO partially contained events. The data (solid circles) 
are seen to fall below the null oscillation expection in every bin of both samples. 

Fig. [7|a shows the three CC event categories studied by MACRO. These include: 
i) "In-Up" events, which are (v^ + V^) reactions occurring inside the detector creating 
muons which exit through the top; ii) "Up-Stop" and "In-Down" events which are 
classified on the basis of topology (timing information not available) and which are 
analyzed together; and Hi) "Up-Through muons" which are initiated by high energy 
(y^ + Vy) interactions below the detector creating muons which traverse the detector 
from bottom to top. Parent neutrino energy distributions for each of the three event 
categories are shown in Fig. |7]b. 
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The In-Up events, and also the Up-Stop plus In-Down events, probe the multi- 
GeV E p regime; the mean E v values for these two MACRO samples are 5 GeV and 
4 GeV respectively. Zenith angle distributions for these partially contained samples 
are shown in Fig. || The event populations are binned in cos6 z from the horizontal 
at cos 9 Z = 0.0, to vertically upward at zenith cosine -1.0. For the 116 events of the 
In-Up sample shown in Fig. |8|a, the data fall below the null oscillation Monte Carlo in 
every bin. (Note that the acceptance for this planar calorimeter is relatively lower for 
horizontal directions.) The ratio of In-Up events observed to the MC prediction, is 
0.57 ± 0.16. Thus the In-Up sample exhibits the large-scale depletion for multi-GeV 
upward going events seen in the Super-K data. The data are seen to distribute in 
accord with the oscillation best fit based upon MACRO Up-Through muons which 
is described in the next Section. A similar trend is observed with the Up-Stop and 
In-Down sample shown in Fig. ^o. Since the latter sample contains roughly equal 
portions of Up-Stop events which are fully oscillating and of In-Down events from 
above horizon which are not oscillating, the amount of depletion relative to null 
oscillation is reduced compared to that of the In-Up sample. 

8 Upgoing Muons in Super-K and MACRO 

There are two event samples which are initiated by below-detector [v^ + z7 M ) inter- 
actions, namely upward stopping muons and upward through-going muons, and they 
represent two different portions of the neutrino spectrum. This can be seen from 
comparison of E v distributions (2) and (3) of Fig. |7|b, which roughly characterize the 
muon samples in Super-K as well as in MACRO. The E v spectrum (2) which produces 
upward stopping muons is distinctly lower, with the bulk of the spectrum lying below 
40 GeV. The different E v regimes give rise to rather different oscillation behavior, as 
can be seen by evaluation of the phase angle in Eq. (4) at our nominal — > v T pa- 
rameter values sin 2 26 = 1.0 and Am 2 = 5 x 10~ 3 eV 2 . At E u = 40 GeV, the vacuum 
oscillation length L Q equals 1.5 Earth diameters. Recall that L is proportional to 
E v and that the oscillation phase is ttL/L . Then for E u much larger than 40 GeV 
as in the case for many through-going muon events, Lq exceeds the Earth's diameter. 
The result is that z/'s initiating through-going muons generally have small oscillation 
phase angles and hence give rise to low oscillation probabilities. On the other hand, 
for up-stopping muons, E v < 40 GeV and neutrino L values are less than L values 
so that sizable phases and large oscillation probabilities frequently occur. 

The expectation that oscillation will occur in relatively different proportions in up- 
stopping versus upward through-going muon samples has been examined by Super-K. 
This is done by measuring the upward-stopping to up-through ratio of muon fluxes. 
In the presence of the two-state mixing inferred from the in-detector samples of Sect. 
6, the muon stop/thru ratio should fall below the MC prediction for null oscillation. 
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Fig. |]a shows the stop/thru ratio plotted versus cos9 z for muons incident from the 
horizontal (cos 6^ = 0.0) to those most vertically upgoing (cos# z = —1.0). The 
observed stop/thru ratio is 0.24 ± 0.02 which is 2.8 standard deviations below the 
null oscillation expectation of 0.37 ± 0.05. 
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Figure 9: a (left): Ratio of upward-stopping to through-going muons versus cos6> z observed in 
Super-K. The data ratio falls below the null oscillation expectation for each of the five below-horizon 
bins, b (right): The upward through-going muon flux versus zenith angle observed in Super-K. 
The shape of the data distribution (open circles) differs from the null oscillation expectation (solid 
histogram) and is better described by two-state oscillations (dashed histogram). 

Additional information can be gleaned from the upward through-going muons 
alone. These events arise from neutrino interactions which have the highest range 
for parent E v . To see how oscillations affect this sample, consider E u at 100 GeV, 
which is the mean energy estimated for Super-K events. (However the distribution 
of parent energies is broad and extends above 1000 GeV.) At our representative 
parameter settings, the vacuum oscillation length is approximately 1.2tt times the 
Earth diameter; consequently the phase angle of the flavor oscillation probability is 
L/1.2 in units of Earth diameter. For horizontal muons the flight paths of parent 
neutrinos are of order 500 km or 0.04 Earth diameters, and so the neutrino phase 
angles will be too small to induce significant oscillation probability. However for 
vertical muons the neutrino paths L become comparable to the Earth's diameter, 
and the phase angles become sufficiently large that rapid oscillation swings ensue. 
( For E v < 100 GeV, oscillations will also occur for muons incident away from the 
vertical. ) 

Available for this conference is an updated through-going muon sample exceeding 
one thousand events from Super-K, the zenith angle distribution of which is shown 
in Fig. ^]b. (This sample is noticeably larger than the one published this spring from 
a 537 day exposure J0J.) Fig. ^b exhibits the trends implied by oscillations for this 
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Figure 10: a (left): The angular distribution of upward through-going muons observed in 
MACRO. The data distribution (solid circles) differs from the null oscillation expectation (shaded 
band) in shape and in rate, b (right): The neutrino oscillation allowed regions obtained by MACRO 
from the upward through-going muons. Confidence level and experimental sensitivity boundaries 
are calculated using the Feldman-Cousins method. 



sample: For bins which are just below the horizon, the data agree with the Monte 
Carlo expectation for null oscillation. For the bins below cos 6^ = —0.3 however, the 
data fall below the null oscillation prediction and this trend continues with muons of 
more vertical inclination. That is, the shape of the zenith angle distribution of these 
upward throughgoing muons, and their overall flux rate as well, deviate significantly 
from null oscillation and agree with expectations from — > u T mixing as inferred 
from in-detector neutrino interactions. 

The same features are seen in the angular distribution of upward through-going 
muons recorded by MACRO [^0J as shown in Fig. |TD|a. The ratio of data to the 
MC prediction for MACRO is 0.74 ± 0.03 ± 0.04 ± 0.12; the last term containing the 
largest uncertainty reflects limited knowledge of the absolute neutrino flux and of 
deep inelastic neutrino cross sections. 

Fig. |il]| b shows the allowed region of the oscillation parameters obtained by MACRO 
based upon the upward through-going muon sample. The allowed region is calcu- 
lated using the Feldman-Cousins method [plfl . For the MACRO data, the minimum 
X 2 point (x 2 of 10.6) is in the unphysical region at sin 2 29 = 1.5 . To clarify the 
situation, the experimental sensitivity at 90% CL is also provided. This is the 90% 
CL contour that would have been obtained had the data coincided with the oscil- 
lation MC expectation at the nearest point inside the allowed region (sin 2 2$ = 1.0; 
X 2 = 12.5). 
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9 Best Fits for sin 2 26 and Am 2 



Available for this Symposium is a new 'all data fit' by Super-K which determines the 
Up — > v T allowed region using a \ 2 summed over all neutrino samples including FC + 
PC events (848 livedays) plus up-throughgoing muons (923 livedays) plus up-stopping 
muons (902 livedays). As shown in Fig. |ll~|a, regions in the sin 2 29, Am 2 plane allowed 
by this fit are the most restrictive ones ever obtained. The \ 2 f° r the oscillation best 
fit is 67.5 for 82 d.o.f., to be compared with 214 for 84 d.o.f. for null oscillation. The 
best fit values are sin 2 29 = 1.0 and Am 2 = 3.5 x 10 -3 eV 2 with flux normalization 
parameter a = +0.06. 
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Figure 11: a (left): Allowed regions obtained by Super-K based upon \ 2 fitting to FC and PC 
single ring events, plus upward stopping muons, plus upward through-going muons. b (right): 
Oscillation parameter allowed regions from Kamiokande (thin- line boundary), Super-K (thick- line 
boundary), MACRO (dashed boundaries), and Soudan 2 (dotted and dot-dashed boundaries). 



In order to gauge the overall consistency of atmospheric neutrino observations, the 
various oscillation-parameter allowed regions obtained by each of the underground 



experiments have been assembled in Fig. |TT|p. Along with the 90% CL region from 
the Super-K all data fit discussed above, we include the final 90% CL region reported 
by Kamiokande (thin-line boundary) J22]. From MACRO, we show the 90% CL 
region and include the experimental sensitivity at 90% CL (dashed boundaries). For 
Soudan 2, we plot the region allowed at 68% and 90% CL. Fig. |Tl"] b shows that, 
for two independent water Cherenkov detectors and for two quite different tracking 
calorimeters, there is a region of oscillation parameter values which is acceptable 
to all experiments. Now "Fools rush in where Angels fear to tread", it is often 
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said. And concerning the significance of Fig. |TT| b the Angels will urge caution, for 
the atmospheric data show neutrino disappearance only - oscillation appearance has 
yet to be shown. Nevertheless, this Fool cannot resist the rush: I propose to you 
that congratulations are in order for the researchers of Kamiokande and of Super- 
K and more generally, for the non-accelerator underground physics community. For 



Fig. lib, Ladies and Gentlemen, is the portrait of a Discovery - the discovery of 



neutrino oscillations with two-state mixing. 



10 Dominant and Sub-Dominant Two-State Mix- 
ing 

Assuming that muon neutrinos oscillate into other flavor (s), with nearly maximal 
mixing and with Am 2 in the range 10~ 2 to 10 -3 eV 2 , it is of interest to consider 
what flavors are involved in the dominant two-state oscillation, and in other possible 
sub-dominant oscillations. That — ► v e could be the dominant mode for disap- 
pearance is ruled out by the CHOOZ reactor experiment. CHOOZ has established a 
limit on V e disappearance ||23|| ; by CPT symmetry, this limit implies that v e neutri- 
nos do not disappear, or at least not in a parameter regime which is relevant to the 
atmospheric flux. 

Since it is generally believed that — > v T oscillation is the dominant mode, it is 
natural to ask: Where are the v T events? In the dominant — > v T scenario, about 
0.9 charged current v T events per kiloton year exposure can be expected to occur in 



an underground detector | 2^| . Then, in exposures reported at this Symposium, we 
would expect Super-K to have recorded ~ 47 FC or PC charged current v T events 
and Soudan 2 to have recorded about 4 events. These events will be up-going but 
otherwise indistinguishable from energetic NC events, and so there is little hope that 
v T reactions can be isolated by the on-going atmospheric neutrino experiments. It 
would be heartening to see a few unambiguous tau-neutrino interactions - even from 
an accelerator experiment! On this, hopes are placed with the v T candidates recorded 



by the DONUT hybrid emulsion experiment at Fermilab [21}]. 

Although dominant — > v T is unlikely to be confirmed anytime soon via v T 
appearance, progress has been made by Super-K towards eliminating the remaining 
competition which is — > v s teriie oscillations. Now sterile neutrinos, by definition, do 
not interact with normal matter via neutral currents, a fact which has consequences 
currently being examined by Super-K. Firstly, it follows that sterile neutrinos cannot 
produce single 7r° events since these are NC reactions: v s N ■/->■ v s N^° . Then, relative 
to the Pn — > v T scenario, Vn — > v s will result in fewer single 7r° events, and the relative 



dearth of these final states will be in up-going directions | 26| . Unfortunately, cross 



sections for these NC reactions have large uncertainties, a situation which hinders the 



isolation of a depletion which is demonstrably significant [27] 
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Figure 12: a (left): Allowed regions for — ► v T and for — > v ster u e two-state mixing scenarios, 
obtained from FC single ring events of Super-K. b (right, top): Zenith angle distribution for selected 
energetic PC events of Super-K (solid circles); histograms show the — > v r and — ► f s solutions 
of Fig. |T^a. c (right, bottom): Up/Down ratio versus log (Am 2 ) for the energetic PC events. 
Curves depict predictions from i/n — > v T (dashed) and Vn — ► v s (dotted) oscillation scenarios. The 
observed ratio ±ler defines the horizontal band allowed by the data. 



More generally, the absence of coupling to the Z° for sterile neutrinos means that 
their effective potential in matter differs from that for neutrinos. The difference 
in respective matter potentials can be written 

r- N 

V,-V s = tV2G f ^, (7) 

where the difference is negative (positive) for neutrinos (antineutrinos) and N n is 
the neutron number density. For high energy z/^'s traversing matter, the existence 
of this potential difference causes — * v s to be suppressed relative to — > v T 



in a way that may be discernible |2^. |29 |. To elucidate the effect, the neutrinos of 



lower energy observed in Super-K can be used to establish a baseline for neutrino 
oscillation behavior. Fig. |T2^ a shows allowed regions in the parameter plane obtained 
from three different oscillation fits to the FC single-ring events. The allowed region 
for Up — > u T oscillations (Fig. |T2]a-top) is as described previously with (1.0, 3.0 x 10~ 3 
eV 2 ) as best values. For — > v s there are two solutions with two allowed regions 



17 



(Fig. |T^a-middle, bottom); these arise from the two possibilities with the sign of 
the mass squared difference which occurs between mass eigenstates involved in the 
mixing. However the allowed regions are found to be very similar in all three cases, 
with fits of comparable quality, consequently — > v s is practically indistinguishable 
from — > v T for the FC events. 

For higher energy neutrinos however this 'degeneracy' can be altered by matter 
traversal, a possibility which can be seen by examining oscillation phenomenology 
appropriate for neutrinos moving through matter of uniform density. Interestingly, 
matter effects for v oscillations can be formulated in a way which is look-alike to 
phenomenology for vacuum oscillations [29|, One feature is that the mixing angle 



for vacuum oscillations sin 20 goes over to sin 20 m for oscillations in matter 

• 2™ sin 2 29 

sin 2 20- 



sin 2 29 + (D- cos 29) 2 ' v 1 

where D is proportional to E v and to the difference in effective potentials in matter 
for the mixing neutrino flavors 

D = 2 -^f where V aP = V a - Vp. (9) 

For z/ M — > u T , V^r = and hence D = 0, consequently matter traversal produces 
no effect on this oscillation. For — *> v s however, and for — > u e as well, V^ s (V^ e ) 
is non-zero, and is in fact sign- dependent since anti-neutrinos are affected differently 
than are neutrinos. Then D is non-zero and acquires a sizable magnitude at high E v . 
Because of its occurrence in the denominator of sin 2 20 m , it acts to suppress — > v s 
at high energies, a suppression which is absent for — > v T . 

To test for occurrence of matter-induced oscillation suppression, SuperKamiokande 
has examined two different high energy neutrino samples. The first sample consists of 
PC events for which the number of photo-electrons from each event exceeds 45,000. 
This is equivalent to requiring that E vis > 5 GeV; the sample thus obtained has 



< E v >~ 25 GeV. The zenith angle distribution for these events is shown in Fig. |12|b. 
The data (solid circles) are binned in cos0 2 from -1.0 to 1.0. Shown superposed are 
the predictions [29| from — > v T and for — > v sy with Am 2 and sin 2 20 set to 
the values inferred from the FC events. A difference between these distributions is 
apparent for cos0 z below -0.2. Neutrinos which initiate events in this region travel 
thousands of kilometers through the Earth, and thus experience matter effects. For 
the Vsteriie oscillation case, matter effects suppress the — > v s oscillation, conse- 
quently fewer z/^'s "disappear". The expectation therefore lies above the curve for 
Up — » u T . Interestingly, it also lies above the data. 

To quantify the difference, an up-down ratio is used. Here, the number of u^s 
which are upward-going (and consequently subject to matter suppression for the v ^ — > 
v s case) is compared to the downward-going flux which, at high energies, is not affected 
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by oscillations: (iV up (cos# < -0 .4) / 'N down (cos 9 > +0A)) DATA = 0.50 ± 0.12 ± 0.01. 
In comparison, the ratio N up /Nd OW n is is 0.94 ±0.04 from the null oscillation MC. As- 
suming sin 2 29 = 1.0, this ratio can be plotted versus Am 2 as in Fig. |T2| c. The region 
allowed by the data corresponds to the horizontal band centered at N up /Nd OW n = 0.50 



with boundaries at ±lo~. Curves obtained from the v.. 



v T and v.. 



v s scenarios 



are also drawn. For v.. 



the predicted curve falls within the band allowed by 



the data for plausible values of Am which include the best fit value obtained with 
the FC events. For — > v s however, the scenario curve lies above the allowed band, 
only "entering" at a Am 2 which is higher than the FC best fit value. 
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Figure 13: a (left, top): Flux of upward through-going muons versus cos 6^ in Super-K. The 
histograms depict — > v T and — > v s solutions, b (left, bottom): Vertical/Horizontal ratio 
versus log(Am 2 ) for muons of 13a, with predictions from — > v T (dashed band) and — > v s 
(dotted band) oscillation scenarios, c (right): Exclusion regions for — > v r and — > v s scenarios 
obtained by fitting to Up/Down and Vertical/Horizontal ratios from energetic PC and through-going 
muon samples, to be compared to the allowed regions of Fig. [l^a. 

A similar pattern is found with zenith angles for upward through-going muons 
in Super-K, a sample for which < E v > ~ 100 GeV. Figure [L3|a shows the cos 6^ 
distribution of that sample, with — > v T and for —>■ v s shown superposed. As 
observed with energetic PCs, matter suppression for z/ M — > v s places the prediction 
(at sin 2 29, Am 2 from the FC events) above the data for angles of incidence corre- 
sponding to large path lengths through the Earth. To quantify the differences in 
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scenarios here, the data distribution is separated into "horizontal" (cos 8 > —0.4) 
and "vertical" (cos 6 < —0.4) events and the vertical to horizontal ratio is calculated: 
(^vertical I ^horizontal) dat a = 0. 77 ± 0.05 ± 0.01 . As shown in Fig. [13|b, there is no value 
of Am 2 for which the v u — > v s scenarios predict N v /Nh in the value range indicated 
by the through-going muon data. 

The difference between data and MC predictions for the three oscillations sce- 
narios can be evaluated using a chi-square function x 2 — Xpc + Xthru-u where X 2 pc 
contains the difference for the up/down ratio and Xthru-u the difference for the ver- 
tical/horizontal ratio. Parameter regions excluded at 90% and 99% CL are then 
deduced from x 2 , corresponding approximately to 2a and 3a exclusion. The exclu- 
sion regions for each of the three oscillation scenarios are displayed as shaded areas in 



Fig. |I~3|c. Comparison of these excluded regions with the oscillation-allowed regions of 

for both Am 2 > and < 0, to be excluded 

/„, it is clear that 



Fig. |12|a, shows large portions of — >• v s 

at the 2a level. While the observations do not as yet rule out u a 



this new approach by Super-K can be steadily strengthened with more exposure. 
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Figure 14: a (left): Triangle graph for displaying possible flavor compositions of the mass eigen- 
state v 3 . b (right) Regions of allowed (U^ 3 , U^ 3 , U^ 3 ) versus Am 2 , obtained by fitting the Super-K 
and CHOOZ data via three-flavor mixing in the approximation of one mass scale dominance pl|. 



11 Subdominant <->• v e in Three-Flavor Mixing 

For a view of possibilities with subdominant oscillations, we turn to investigations of 
three-flavor mixing. A number of approaches have been discussed in the literature 
PH. Here we review a scenario for subdominant — > u e which emerges directly from 
the approximation of one mass scale dominance. In this approximation it is assumed 
that one of the mass eigenstates - v% let us say - is more massive than the other two, 
and that the lighter eigenstates are nearly mass-degenerate. As a result there 
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are two mass-squares differences Am 2 



m 



2 2 1 and Sm 2 



\m% 



m; 



having 



different magnitudes. The larger Am 2 can be identified with the dominant two- 
state mixing observed with atmospheric oscillations, whereas the much smaller 5m 2 
relates to oscillations not observable with atmospheric z/'s but presumably relevant 
to solar z/'s. Up to terms of order (5m 2 ) j (Am 2 ) the parameter space for atmospheric 
neutrinos is spanned by (Am 2 , U 2 3 , U 2 3 , U 2 3 ). The amplitudes satisfy the unitary 
constraint U 2 3 + U 2 3 + U 2 3 = 1. 

For vacuum oscillations, it follows that the oscillation probability for transitions 
between flavors a and /3 is 



= 4 U; 



a3 ■ U p3 ■ sin' 



1.27 Am 2 ■ L 

R 



(10) 

As suggested by the amplitudes in this expression, the flavor composition of the 
massive u 3 eigenstate is the central issue: 



^3 = U e3 ■ v e + ■ + U t3 ■ v T . (11) 

For fixed Am 2 , the v 3 composition is conveniently depicted using the equilateral trian- 
gle construction shown in Fig. Ilia for which the unitarity constraint is automatically 



incorporated. 

Fogli, Lisi, Marrone, and Scioscia have compared predictions for specific choices 
of (Am 2 ,U 2 3 ,U 2 3 ,U 2 3 ) with Super-K zenith angle distributions; the constraint on 
v e — * from the CHOOZ limit has also been included. They find that two-flavor 
oscillations with maximal <-> v T mixing works rather well [ |31|j : 

sin 2 26^ = 4:U 2 3 - U 2 3 ~ 1.0 with U 2 3 ~ U 2 3 ~ 1/2 and U 2 3 ~ 0. (12) 

A small admixture of U 2 3 is however allowed by the fits to data as is shown graphically 
in Fig. |j"4|b. Here, for a relevant selection of Am 2 values, the domain of U 2 3 values 
allowed by Super-K data at 90 and 99% CL comprise the shaded areas in the triangle 
graphs of the left-most column. Elimination of regions excluded by CHOOZ (see 
center-column triangle graphs) leaves the diminished but still existent allowed regions 



shown in the right-most column of Fig. 14b. From the height of the various allowed 



regions, it is concluded that U 2 3 ~ 0.05. The expression for the (z/ M — > u e ) vacuum 
oscillation probability follows immediately from Eq. (|T0|) with a, (3 assigned to /x, e 
respectively. We infer from this formula that P mc (z/ /i — > u e ) can be as large as 0.10. 



12 P(V[i —> v e ) Amplification via Matter Resonances 

A <-> v e oscillation of strength as indicated above will be hard to discern within 
the atmospheric flux, however we may get some help, as a result of amplification by 
matter resonances in the Earth. Two kinds of resonance effects are possible. 
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Figure 15: a (left): Values of sin 2 29 and E v for which to v e mixing of atmospheric neutrinos 
could be enhanced by MSW resonance in the Earth's mantle or core P^| . b (right): Possible "bump" 
in up/down asymmetry of v e flavor neutrinos which could arise with MSW resonant enhancement 
in the Earth's mantle |B3l. 



One arises with the MSW resonance wherein the D of Eq. (|9]) containing V eil ap- 
proximates cos 2 26 with the consequence that sin 2 26 m ~ 1.0. Depending upon the 
particular values of the mixing parameters, MSW enhancement can take place either 
in the terrestrial mantle or core for the E v intervals depicted in Fig. |T5| a [[32]. An 



MSW resonance could result in a bump in the upward-going v e flux at the resonance 
energy. This possibility has been examined by J. Pantaleone who proposes that the 
v e up-down asymmetry could be a useful discriminant, with possible outcome as il- 



lustrated in Fig. |i5|b [33]. Of interest to long baseline experiments, e.g. K2K and 



MINOS, is the observation that MSW enhancement can also take place in the Earth's 



crust 34 



A second and different resonance-like enhancement can take place for atmospheric 
neutrinos which cross the Earth's core. For such neutrinos, having paths that cross the 
mantle, the core, and again the mantle, a complex constructive interference among the 
oscillation amplitudes arising in regions of different density is possible. The algebraic 



delineation of this effect has been presented by M.V. Chizhov and S.T. Petcov [35]. 
Their work has yielded striking depictions of the transition probabilities as shown in 
Fig. |16[ (Alternative formulations and interpretations have been presented; see Refs. 

Fig. [I~6| a shows v e <-> oscillation probability as a function of nadir angle (neutri- 
nos at vertically upward incidence being at 0°) and of E u /Am 2 in units of (MeV/eV 2 ). 
Here the smaller peak at E^/Am 2 = 2.5 x 10 6 is the MSW resonance in the Earth's 
mantle. The distinctly larger structure arises with the mantle-core-mantle trajecto- 
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Figure 16: a (left): The probability Y{y e ^ — > ^ M;T ( e )) for atmospheric neutrinos with mantle-core- 
mantle trajectories through the Earth. Such neutrinos may undergo a resonance-like enhancement 
(the large peak) which is different from an MSW enhancement, b (right): Regions of resonance-like 
enhancement for Earth core-crossing neutrinos |35|] . 

ries. The structures shown in Fig. [TE| a are to be found on the "probability island" at 
lowest sin 2 29 in Fig. |16|b. Other resonance structures exist at higher mixing angle 



values as shown. 

In the on-going underground experiments the effective integrations over E u , cos 9 Z , 
and detector resolution effects, which necessarily occur with the accumulation of data 
events, likely assure that matter resonance effects will be difficult to observe. Be that 
as it may, these are intriguing phenomena, underwritten by phenomenology which is 
rich and well-grounded. Their elucidation poses an interesting challenge for a neutrino 



factory to be built at a muon collider 37 



13 Atmospheric Fluxes; Concluding Remarks 

Invaluable to all oscillation analyses with atmospheric neutrinos are developments 
which yield improved knowledge in rates and shapes of atmospheric flux spectra. 
There have been several such developments of recent, which we briefly remark upon 
here. Firstly there is the observation by Super-K of the east-west anisotropy in 
horizontal neutrino fluxes at the Kamioka site |^9f| . The and v e fluxes from the east 
are found to be depleted, as expected due to geomagnetic cutoff of charged cosmic 
ray primaries and as calculated in the one-dimensional models of the atmospheric 



cascade HOI. Just arrived "on the scene", are three-dimensional atmospheric flux 



calculations which have been prepared independently by two groups [[41]]. Their arrival 
is timely indeed, since 3-D calculations are the natural framework in which to utilize 
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the abundant data becoming available from balloon-borne spectrometer experiments. 
These experiments measure the primary cosmic ray flux and sample the secondary 



muon fluxes at a variety of depths in the atmosphere [42 



Concerning flux-related measurements, there are "swords-in-the-stone" aplenty to 
tantalize the brave-of-heart. For example, no experiment to date has separated and 
compared the atmospheric anti-neutrino fluxes to the neutrino fluxes. It is within the 
capability of the underground experiments to distinguish z7 M from interactions Jl3| . 
and the v, V reaction cross sections are known. Nevertheless, the assertion by all of the 
flux calculations that the ratio of ^(i 7 ) to is very nearly 1:1 for either atmospheric 
neutrino flavor, remains untested. To this end, a V^jv^ ratio of ratios measurement 



would be interesting |[44|| . As a second example, we note the absence of measurements 
which examine variations in the atmospheric neutrino fluxes predicted to occur as a 
function of the solar cycle. The variations are substantial in the sub-GeV portion 
of E v spectra and should be more pronounced at northern geomagnetic latitudes 



4q| . Such a measurement of course places a premium on continuous exposures which 
extend to a decade or longer; however 1MB, Kamiokande, and Soudan 2 have shown 
solar-cycle-duration exposures to be attainable. 

To conclude: A bonanza in neutrino oscillations research is in progress, driven to 
fever pitch by recent experimental observations with atmospheric neutrinos. Quite 
possibly, the finest nuggets are still in the ground. Fortunately the atmospheric beam 
is always on and beamline access is free; however the detectors required for future 
progress will not materialize cheaply In any case, the aura of adventure and 

discovery which now pervades the atmospheric neutrino beamline will remain for 
some time. The Dreamers and the Restless will come to try their luck, and among 
them - appearance probability of 1.0 - will be participants from this Symposium. 



It is a pleasure to thank the organizers of the Symposium, and especially Helen 
Quinn and John Jaros, for the opportunity to give this Talk. I am greatly indebted 
to Takaaki Kajita, Serguey Petcov, Ed Kearns, John Learned, Francesco Ronga, 
Maurizio Spurio, Tomas Kafka, Jack Schneps, Maury Goodman, and Sandip Pakvasa 
for communications and discussion relating to physics with atmospheric neutrinos. 
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Discussion 



B.F.L. Ward (University of Tennessee): How do we combine the results for Am 2 
from SuperKamiokande and for Soudan 2, for example 3.5 x 10~ 3 eV 2 and 8 x 10~ 3 
eV 2 ? 

Reply: I would not recommend doing that. The Soudan 2 and MACRO measure- 
ments are interesting as checks, with completely different technique and systematics, 
on the Super-K result. But since the three determinations are in agreement and 
since the Super-K measurement is the one with predominant statistical weight, the 
Super-K value is the one to be used. 
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Peter Rosen (DOE): The evidence for oscillations from atmospheric neutrinos is 
certainly impressive and the community of non-accelerator physicists is certainly to 
be congratulated. To what extent can you rule out alternative explanations? For 
example, Vernon Barger, Sandip Pakvasa et al. have shown that the data can be 
fitted by a neutrino decay scenario. 

Reply: The work to which you refer |I7) makes a good case for neutrino decay being 
viable as an alternative to neutrino oscillations. Also, there are other explanations, 
e.g. flavor changing neutrino interactions, which are not ruled out [E3 . 



S. Ragazzi (University of Milano): What do you expect to learn from the com- 
parison of $(z7) and $(z^) ? 

Reply: The atmospheric neutrino flux calculations predict anti-neutrino fluxes to 
be nearly the same as neutrino fluxes; this should be tested. Although no difference 
is to be expected from the viewpoint of conventional flavor oscillations, I note that 
neutrino into anti-neutrino oscillation schemes have good lineage, originating with 
Pontecorvo's proposal of 1957 [49 . 



Jasper Kirkby (CERN): Have you looked in your data for a signal of 7's produced 
by solar cosmic rays? These are produced by events lasting a few days and created 
by energetic coronal mass ejections from the sun. They produce particles with peak 
energies of about 100 MeV, and occasionally up to about 1 GeV. During the most 
energetic events a large ionization — equivalent to 20-30% of the total annual galactic 
cosmic ray flux — is dumped into the Earth's atmosphere. These events occur near 
solar maximum, which we are entering now. I would guess Super-K may be able 
to detect i/'s in-time with these events. Perhaps they could even contribute to the 
distortion of the solar neutrino energy spectrum we saw in the previous talk in the 
hep v energy region. 

Chang-Kee Jung (SUNY, Stony Brook): Super Kamiokande has examined data 
for solar activity dependence for long-term periods. But we have not done so for 
specific short period dependence on solar flares. 



28 



